The binding dynamics of tropomyosin on actin 
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We discuss a theoretical model for the cooperative binding dynamics of tropomyosin to actin 
filaments. Tropomyosin binds to actin by occupying seven consecutive monomers. The model 
includes a strong attraction between attached tropomyosin molecules. We start with an empty 
lattice and show that the binding goes through several stages. The first stage represents fast 
initial binding and leaves many small vacancies between blocks of bound molecules. In the second 
stage the vacancies annihilate slowly as tropomyosin molecules detach and re-attach. Finally the 
system approaches equilibrium. Using a grain-growth model and a diffusion-coagulation model 
we give analytical approximations for the vacancy density in all regimes. 



INTRODUCTION 

Tropomyosin (TM) is the protein that plays t he key 
role in the regulation of muscle contraction (Lchrcr 
and Gccvcs, 1998| ) . It comes in a rod-like form and 
binds to the groove each side of double-helical actin 
filaments (Fig. [|). Actin filaments give the cytoskele- 
ton its mechanical stability and serve as tracks for mo- 
tor proteins from the myosin family, especially in mus- 
cle cells. In mus cle each tropomyosin molecule covers 7 
actin monomers ( Hitchcock-deGrcgori and Varnell, 1990 



McLachlan and Stewart, 1976). Other tropomyosins 



covering 6 and 5 monomers exist as well (reviewed by 
Hitchcock-deGrcgori (199~4| )), but for specificity we will 
focus on TM covering 7 actin subunits. Together with the 
associated protein troponin, TM on actin can switch be- 



twee n two laterally shifte d conformations ( Lehman et aL 
200(j; Vibert et aL, 1997 ). The equilibrium between these 



two conformations is strongly influenced by the concen- 
tration of Ca 2+ ions. As one confor mation obstructs the 
stron g binding of myosin to a ctin ( |Geeves and Lchrcr 



199|; Lehrer and Geeves, 1998 ), this provides the mech- 



anism of calcium-mediated activation of myosin in skele- 
tal muscle cells. The regulation of myosin activity has 
been widely studied using Actin-tropomyosin fi laments 
assembled in vitro (Fraser and Marston, 1995; Gordon 
Other functions of TM might include the 



et a 



1997) 



stabiliza tion of actin aga i nst fragmentation and actin as- 
sembly ( Lazaridcs, 1976 Wcigt et aL, 1991 ). 

The binding of t ropomyosin to ac t in has been st ud- 
ied in equilibrium ( Hill et aL, 1992 ; Wegner, 1979) as 



well as dyn amically ( Wegner and Ruhnau, 1988| ; Weigt 



et al, 1991) and has been found to be highly coopera- 



tive ( |Hill et aL, 199$ |Wegner, 1979J ). But the way that 
tropomyosin binds to actin is not simple. As each TM 
molecule occupies 7 actin monomers, an obvious problem 
is that sometimes a gap of 1-6 monomers, too short for 
another TM molecule to be inserted, can remain on actin. 
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FIG. 1: Tropomyosin (TM) binds on actin poly- 
mers on each side of the double helix. One 
TM molecule thereby occupies 7 actin monomers. 



It has been suggested (Weigt et aL, 1991) that these gaps 
could play an important role by regulating the binding 
of a-actinin and the fragmentation kinetics of actin. 

In this Article, we study a theoretical model for the dy- 
namics of binding of TM to actin based on the reaction 
kinetics sh own in Fig. |2], first proposed by Wegner and 
coworkers (Wegner, 1979; Weigt et aL, 1991). The equi- 



librium distribution of bound molecules and vacancies 
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has been known exactly for a long time ( Wegner, 197£ ) 
and repre sents a special case of a more general model 
solved by McGhee and von Hippel (1974 ). However, as 
we have pre viously shown in the context of kinesin on mi- 
crotubules ( Vilfan et aL, 2001 c ,^), the time necessary to 
reach the equilibrium state can be very long. There are 
many similarities between the binding of dimeric kinesin 
on microtubules and TM on actin. In both situations 
there is an initial period in which all gaps that are wide 
enough to accept a new molecule get filled. The dynamics 
within this period has been directly experimentally ob - 
served using fluorescent techniques ( Weigt et aL, 1991 ). 
After the initial period, only vacancies stretching over 1 
to 6 actin monomers, too small for another TM molecule 
to fit in, remain. They get healed out in a much slower 
annealing process in which bound molecules detach and 
others re-attach at other places. This process takes place 
on a much slower time scale. The annealing process is 
finished when the system reaches the equilibrium state. 
Despite the similarities with the binding of kinesin onto 
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FIG. 2: Reaction scheme for all binding and unbinding 

processes of TM on actin. The binding (unbinding) rate 

(k_), at the "— " side of an- 
M„ (U M\ at the «,« gide 



for an isolated site is k+c 
other bound TM molecule kf c (k M ) 
k+c (fc^), and between two bound TM molecules k\ 



denotes the concentration of TM in solution. 



microtubules, there are also important differences: 

1 . The TM molecules occupy 7 instead of two binding 
sites. 

2. Kinesin dimers can bind with one head as well as 
with two, while TM must either bind all sites or be 
fully detached. 

3. The interaction between bound tropomyosin 
molecules is much stronger. 

In our theoretical analysis we make reasonable simplifi- 
cations justified by experimental observations and then 
use results from diffusion-annihilation models to find an 
analytical solution for the dynamics in our model. The 
central result is the prediction of the number of defects 
(gaps between contiguous blocks of tropomyosin) as a 
function of time for a given solution concentration. 

The fact that a tropomyosin molecule occupies 7 bind- 
ing sites makes it especially interesting from the theo- 
retical point of view. The size 7 lies between the simple 
monomeric or dimeric and continuous models. We have 
shown previously that the relaxational behavior of the 
dimer adsorption-desorption model can be explained by 
a mapp ing onto the A + A — » diffusion-annihilation 
model ( Vilfan et ai, 2001 q ). The gap density decays 
with a power-law oc t~ l l' 1 before it approaches the final 
equilibrium value. This behavior remains qualitatively 
valid even if there is an attractive interaction between 



adsorbed dimers. If one introduces a strong diffusion 
of adsorbed dimers the power-law decay is preceded by 



a mean-field r egime in which ng oc t (Privman and 



Nielaba, 1992). On the other hand, non- interacting k- 



mers (k > 3) are always well described with the mean 
field kinetics of the model kA — > which predicts a ga p 
density decay oc x ) (Nielaba and Privman, 1992). 

In the continuous limit, also called continuous sequential 
adsorption (CSA) or "car parking problem" the length 
of unoccupied space decays as oc 1/ log t in the mean- 
field regime (|Jin et ai, 1994; Krapivsky and Ben-Nairn 



1994). We will show that our model (k = 7 with strong 



attractive interaction) differs from all models mentioned 
above and that it maps onto & A + A ^ %A model, which 
includes particle annihilation as well as coagulation. 



DEFINITION OF THE MODEL 

Our model is based on the reaction scheme shown in 
Fig. |^ and is defined as follows. We assume that there is 
no interaction between bound TM molecules on the op- 
posite sides of the actin polymer and therefore describe 
actin as two independent one-dimensional binding lat- 
tices. We further assume that TM molecules can bind 
to actin only with their whole length, thereby occupy- 
ing 7 binding actin monomers. The rate at which a TM 
molecule binds to a group of seven free binding sites on 
actin is k+c where c denotes the TM concentration. The 
rate for the reverse reaction, i.e., the detachment of a TM 
molecule without neighbors, is k_. The rate at which 
new TM molecules bind beside already bound ones (we 
call such sites single-contiguous binding sites) will be de- 
noted as k+c (the rate at which the new TM binds on 
the plus side of a previously occupied block) and k¥ c (on 
the minus side). The reverse reaction rates are k p and 
k_. It is reasonable to assume that the binding constants 
are only influenced by nearest neighbors. Finally the at- 
tachment rate in a gap of exactly 7 lattice sites will be 
k+c and the detachment rate in the middle of a contigu- 
ous block kR. The principle of detailed balance states 
that the ratio between the binding and unbinding rate 
depends only on the free energy difference and therefore 
has the same value regardless whether a molecule binds 
on the plus or minus side of an occupied block 
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Here J denotes the coupling energy between two TM 
molecules on actin and 7 = exp(J/fcsT) the coopera- 
tivity coefficient (Hill, 1985). The coupling energy from 
both ends is additive, therefore 



kj 
k D 



K 1 



k" 2J 



(2) 



3 



In the following we will assume a strong interaction 
between bound molecules, 7 » 1. The assumption is 
well justified since the cooperativity coefficient has been 
determined to be between 7 = 600 a nd 7 = 1000 b y 
Wegner (19791 ) and around 7 = 100 by |Hill et al. (19921 ). 
The binding affinity of an isolated TM mol ecule has been 
found to be very low (Weigt et at, 1991) and we may 
assume 



K N c « 1 



(3) 



When discussing the dynamics we will neglect all pro- 
cesses where a TM molecule detaches in the middle of a 
contiguous block, setting kP w 0. 



RESULTS 

The equilibrium solution of the problem on an infinite 
lattice has been known exactly for a lo ng time ( McGhcej 
and von Hippel, 1974 ; Wegner, 1979 ). The number of 
attached TM molecules per lattice site tt-tm (which can 
assume values between and 1/7) is given by the follow- 
ing implicit equation 



K N c = 



with 



2(7-l)(l-7n TM ) 
1 - 7n TM V ( 2 7 - - 7n T ]vr) + "tm - R 

1 2(1-7*™) y (4) 

1 -8n TM + RJ V ' 



R = vQ - 8n T ivi) 2 + 477i T m(1 - 7n TM ) 



The average gap size was determined as 

__ 2( 7 -l)(l-7n TM ) 
9 -(l-6n TM )+R 



(5) 



The fraction of unoccupied lattice sites is given as rig q = 
1 — 7riTM and the gap density as n^ q = riQ^/g. In the 
limit of very high TM concentrations and nearly full cov- 
erage, nxM — * 1/7, the average gap size becomes g ~ 1. 

However, the exact equilibrium configuration is of lit- 
tle value without the knowledge of the time the system 
needs to equilibrate (relaxation time). Figure || shows the 
time dependent filament coverage and gap concentration 
as obtained from a dynamic Monte-Carlo simulation for 
different TM concentrations with realistic values for the 
kinetic constants. As the simulation results show, the dy- 
namics shows an interesting two-stage behavior at high 
TM concentrations. In the first stage, TM molecules bind 
to actin wherever there is enough unoccupied space. The 
first stage is followed by a plateau and the second stage 
in which reordering of bound TM molecules takes place. 
In this second stage gaps that inevitably remain after the 
first stage are healed in an annealing process. The sec- 
ond stage is finally followed by an exponential relaxation 



into the equilibrium state given by Eq. |j. The second 
stage is much slower than the first one. We therefore 
conclude that in many experimental situations the rele- 
vant properties of the system are actually determined by 
the dynamic behavior of the model and not its equilib- 
rium state. 



Initial binding 

In this Section we describe the initial dynamics of TM 
starting with an empty actin filament. Since a single 
TM molecule only weakly binds to actin (Eq. ||) nucle- 
ation is needed to initiate binding. F or its description we 
use the approach first developed by Kolmogorov (1937 ) 
in his 'g rain growth' mo del (reviewed by Evans (1993 ); 
see also (Privman, 1997)). Kohnogorov's model assumes 
that nucleation is sufficiently slower than grain growth, 
meaning that grains can grow to large sizes before the 
growth process is stopped by another grain. It therefore 
neglects fluctuations in the growth process and assumes 
a constant growth velocity. The nucleation, on the other 
hand, is described with its full stochasticity. 

Where the blocks growing from two neighboring nuclei 
meet, there is a probability of ^ that they will match 
exactly without leaving a gap (Fig. |J). In all other cases, 
occurring with a total probability of | a gap of 1-6 sites 
will remain. The number of such gaps should therefore 
equal | of the total number of nuclei. All gap sizes be- 



tween 1 and 6 sites occur with equal probabilities and 
the average gap width is 3.5 sites. 

We calculate the number of gaps after the initial bind- 
ing process in the following way. For simplicity we as- 
sume that the actin concentration is low enough that the 
TM concentration in the solution does not change sig- 
nificantly in time (a generalization for a time-dependent 
solution concentration is discussed later). Then the nu- 
cleation rate r n and the growth velocity v are constant 
in time. The probability that a nucleus will be formed 
on site i at time t equals the nucleation rate at that time 
multiplied by the probability that the site has not been 
occupied previously. The latter equals to the probabil- 
ity that there is no nucleus in a triangle of height t and 
base vt in the position-time diagram (akin to the light 
cone frequently used to illustrate independent events in 
the special theory of relativity) (Fig. ||). To improve 
the accuracy at relatively high nucleation rates, we can 
take into account that the nucleus already has a width of 
e?o = 14 sites and the triangle then becomes a trapezium 
with side lengths 2d + vt and 2d and height t. 

The probability that no nucleation event takes place 
within the trapezium is given by the zeroth term of a 
Poisson distribution with the expectation value N — 
r n (vt/2 + 2do)t, namely Pq = exp(-N). Then the nucle- 
ation rate at time t reads 



f n (t)=r n P (t) = r n e- r ^ +2do > 



(6) 



The density of gaps per lattice site in the state after 
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FIG. 3: Gap density (nc(i), dashed line) and fraction of unoccupied binding sites (no(t) = 1 
solid line) as a function of time, obtained from a stochastic simulation. Parameters: 



7titm(*) = n G (t)g(t), 



k 



N 



fef = 50^M _1 5 -1 , kf = 3.0fM.- 1 s- 1 , k N = 100s- 1 and 7 = 1000. The four graphs show data for differ- 
ent TM concentrations: 0.1, 0.25, 1 and lOuM. With these values the effective hopping rate is rh op = 0.28 a -1 . 
The thin lines show the predictions of Eq. [TJ] for the initial phase, Eq. [7] for the plateau, Eq. for the algebraic 
regime and eventually the exact result of Eq. |] for the equilibrium. The agreement between the theoretical calcu- 
lations and the simulation becomes very good for concentrations c > 1 fjM.. Simulation data were obtained on a 
lattice large enough that the results correspond to an infinite system (2 17 sites and periodic boundary conditions). 



initial attachment is then given as 

n 6 



7 Jo 

6 [tt 



dtr n {t) 



r) e 2d ° n Erfc (d Q y/Zq) with r) 



7 V 2 v • / v 

and is plotted in Fig. ^. Its asymptotic limits are 
6 



(7) 



for 
for 



77 — > 00 

ri -h. . (8) 



Once we have expressed the gap concentration with 
the nucleation and growth rate, we need to express these 
with the original model parameters. The total growth 
velocity v (sum of growth velocities at the + and the — 
end, expressed in lattice sites per time unit) reads 

v = 7 [{kl + fcf )c - (k p + k M )} a 7{kl + k?)c . (9) 



In this equation we neglected the effect that uncontigu- 
ously bound molecules (which detach after a short time) 
would have on the growth by blocking the binding sites. 
This approximation is justified as long as K N c <C 1. 

If Kc 3> 1, two molecules bound beside each other 
already form a stable nucleus, while a single one does 
not (Eq. [|). The nucleation rate r„ can be obtained as 
the attachment rate of the first molecule multiplied by 
the probability that a second one binds to its side before 
the first one detaches and reads 



r — k N 



(feg + fcjV 

(fc^ + fcf )c + k 



Together with Eq. [| we obtain 



V 



7({kP + kl I )c + k N ) 



(10) 



(11) 



In the limit of fast equilibration of isolated TM molecules 
(i.e., if the detachment of an isolated molecule is much 
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FIG. 4: Nucleation and growth of tropomyosin (TM) 
blocks on actin. Each time the blocks from two differ- 
ent nuclei meet, a gap remains in 6 out of 7 cases. A 
lattice site at time t can be considered as empty if nu- 
cleation has not taken place in a trapezium of height t, 
base vt + 2do and top 2do where v is the sum of growth 
velocities at both ends and do the size of a nucleus. 



show good agreement with the simulation value as long 
as the assumption Kc 3> 1 is fulfilled. We can go even 
further in our analysis and give an approximation for the 
time dependent vacancy concentration. We can estimate 
the fraction of empty binding sites at time t as the prob- 
ability that no nucleus has reached that point plus the 
number of empty sites that remain in gaps between nuclei 
(on average 3 per nucleus). Therefore, we obtain 

n (t) = 1 - 7n T M(i) = 




(13) 

The comparison between the prediction of Eq. [l3| and the 
simulation result can be seen in Fig. |. The agreement 
becomes very good at concentrations c > 1 fiM. 



Annealing of gaps 
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FIG. 5: The gap density (number of gaps per lat- 
tice site) after the initial binding phase, tlq as a func- 
tion of the parameter 77 = r n /v, as given by Eq. (]?]). 

faster than the attachment of another one to its side) it 
simplifies to 

r) = K N c/7 for k p + k M < (k P + k?)c < k N . 

(12) 

The values of the plateau gap concentration obtained 
from Eq. [?] are shown in Fig. |^ (short dashed line). They 



By now we have calculated the gap concentration after 
the initial binding phase. What follows is a process on 
a much longer time-scale in which molecules at edges of 
the gaps can detach and re-attach. If they re-attach on 
the other side of a gap, the gap makes a move of 7 sites in 
one direction (Fig. ||a) . The rate of such diffusive steps 
to each direction is given as the detachment rate of a TM 
molecule on e.g. the minus side of a block multiplied by 
the probability that a molecule re-attaches on the plus 
side of the other block before another one re-attaches to 
the position where the first one has detached: 

/ m ^+ 1 k p k^ / 1 1 \ 

rhop = - k^ + kf = Kk^ + W = \k F+ k M ) ' 

(14) 

If two gaps, each containing 1-6 sites with equal prob- 
abilities, come together, they can either join to a single 
gap (Fig. ||b) or, if their total size exactly fits one TM 
molecule, annihilate (Fig. ^c). If the original gap sizes 
are g\ and 92, then the joined gap has the size (g\ + 172) 
mod 7. All possible combinations of 171 and gi are shown 
in Table |[ This table shows that if the probabilities of 
both gap sizes g\ and 32 are equally distributed between 
1 and 6, the same will hold for the size of the resulting 
gap g.j. In 1 out of 6 cases the two gaps will annihi- 
late, while they will form a new gap with the same size 
distribution in 5 out of 6 cases. Thus we can represent 
gaps as particles A, hopping randomly along the lattice 
and joining or annihilating when two of them meet. We 
can therefore map our model to a diffusion-annihilation 
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FIG. 6: Examples of effective reactions that occur af- 
ter the detachment and attachment of a tropomyosin 
molecule: Diffusive step (A); Pair coagulation A + 
A -> A (B); Pair annihilation A + A -> (C). 



92 





1 2 3 4 5 6 


1 


2 3 4 5 6 - 


2 


3 4 5 6 -1 


9i 3 


4 5 6 - 1 2 


4 


5 6-123 


5 


6-1234 


6 


- 1 2 3 4 5 



TABLE I: Size of the joint gap after coagulation 
of gaps with sizes g± and gi. If the initial gap 
sizes <?i and g2 are randomly distributed with val- 
ues between 1 and 6, the probability for annihila- 
tion is 1/6. The probability that the joint gap will 
have a certain size between 1 and 6 is always 5/36. 



model consisting of following reactions 



A- 
A 



A 
A 



A 




probability 5/6 
probability 1/6 



According to Lee (1994 ) all diffusion- annihilation mod- 
els of the type 2 A — > I A (or generally kA — > I A) belong to 
the same universality class. In the asymptotic limit, the 
"particle" concentration can be related to that of the ex- 
actl y solvable A+A — > model (Lushnikov, 1987; Forney 
and McConnell, 1983 ) and reads 



n(t) 



1 



I y/8irr hop t 



(15) 



The 



In our case we have to set 1 = 5/6 and r^op = 7 r hop . 
second relation results from the fact that in each diffusive 



step a gap jums over seven sites. We finally obtain 
6 



n G {t) 



49 v /27rr hop i 



and n (t) = 3.5 nc{t) ■ (16) 



Interestingly, the asymptotic particle concentration is in- 
dependent of its initial value. The gap concentration at 
long times therefore does not depend on the intermediate 
gap concentration n° G . As our equation reveals, it is even 
independent of the solution concentration c. 

Of course, the mapping to the A + A — > IA model 
is only valid as long as the concentration of particles A 
is well above its equilibrium value. When these come 
closer, events of pair creation like A — > A + A become 
relevant. One can therefore estimate the equilibration 
time as the time when the vacancy concentration deter- 
mined by Eq. [l6| becomes equal to the equilibrium con- 
centration from Eq. ^. This consideration is only valid 
if the filaments are long enough that even in equilibrium 
a considerable number of gaps remains. Otherwise the 
equilibration time can be estimated as the time in which 
the gap concentration drops below one per filament. 

Figure [5] also shows the comparison between the re- 
sult of Eq. [l6] and the simulation result. The results 
agree well for high concentrations where the intermedi- 
ate plateau and the final equilibrium are far apart. At 
very high concentrations (not shown, as they would be 
beyond experimental relevance), there would be an ad- 
ditional power-law regime between the t~ x / 2 law shown 
here and the final equilibrium. In that regime processes 
of the type A — > A + A become relevant, but not yet 
processes of the type — > A + A. The model then be- 
comes equivalent to the non-interacting 7-mer disposition 
model in which the vacanc y density decays according t o 
the mean- field law cx i -1 / 6 (Nielaba and Privman, 1992). 
But with experimentally relevant parameters one can see 
only a slight remnant of this regime. 



Finite systems 

For now we have assumed infinitely long actin fila- 
ments, neglecting all boundary effects. However, in many 
cases the finite length of actin filaments plays an impor- 
tant role. For example, the length of an actin filament in 
skeletal muscle is just about 1.1 ^m (200 monomers). 

There are different possible scenarios how TM 
molecules should behave at the end of an actin fila- 
ment (Fig. |?]). A hard boundary would mean that a TM 
molecule cannot bind with any segment overhanging the 
end of the actin filament. A soft boundary, on the other 
hand, would mean that binding of a TM molecule is pos- 
sible, albeit with a lower affinity, even if less than seven 
actin sites are free at the end of a filament. The binding 
constant of a TM molecule partially overhanging the end 
of the actin filament can be estimated in the following 
way. If we neglect the entropy gain resulting from the 
flexibility of the overhanging TM end the free energy dif- 
ference between the molecule bound wholly on actin and 
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"The two columns show two possible inerpretations of the exper- 
iment, one with a very low K N and one with a very high K N . 

''Fields for the plateau gap density and equilibration time are left 
empty as no kinetic data is available. 

TABLE II: Model parameters used in our simulation (Fig. |J) and their experimental values. The upper block 
shows the input parameters, the second block the quantities that do not depend on TM concentration, i.e., the 
hopping rate rh op and the time to.oi ( ^o.ooi) when the gap density reaches 0.01 (0.001) per lattice site. The 
lower three blocks show the concentration dependent quantities for three different TM concentrations. These in- 
clude the equilibrium gap density n® q , the average gap size g, the plateau gap density m? G , and the equilibra- 
tion time t Eq . Fields for the gap density are left empty where the equilibrium state is only sparsely covered. 



one with m overhanging segments is 



<5G m — SGq — — Jo 



(17) 



where Jo is the binding energy (including hydrophobic 
contributions) of a TM molecule, which we estimate as 
Jo k 20 ksT. The binding constant for a TM molecule 
with 7 — to bound and to overhanging segments is then 



K„ 



ifexp 



toJ 



7k B Tj 



(18) 



giving the values K x « 0.057 K, K 2 « 0.0033 if, etc. 
With the values we use (K = 10 /iM" 1 and c < 10 //M) it 
turns out that only TM molecules overhanging with one 
or at most two segments can bind. We therefore conclude 



that a hard boundary is a good approximation of the real 
situation, even if it might not be exact. 

The process of relaxation towards equilibrium on fi- 
nite filaments differs from that on infinite ones in several 
aspects. These include the annihilation of gaps when 
they reach the boundary, the creation of new gaps at the 
boundary and the adjustment of the whole TM block to 
reduce the gaps at the boundary. A discussion of finite- 
size effects in diffusion-annihilation models can be found 



in flKrebs et al., 1995Q . 

Figure |§| shows the gap concentration as a function 
of time for different filament lengths, compared with in- 
finitely long filaments. The gap concentration does not 
differ significantly from infinite filaments before it reaches 
a value of about 1 gap per filament length. After that, 
the finite-size effects can lower (due to faster gap anni- 
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FIG. 7: Two scenarios for the boundary condition at the 
ends of an actin filament: a) a hard boundary implies 
that no TM molecule can bind beyond the end of the 
actin filament; b) a soft boundary allows TM molecules 
to bind to the actin end with few segments overhanging. 
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FIG. 8: Average values of the gap density (nc(i), dashed 
lines) and fraction of unoccupied binding sites (no(t) = 
1 — 7n,TM(t), solid lines) as a function of time for fil- 
aments of finite length (200, 1000 and 5000 subunits), 
compared with results on "infinite" filaments, as obtained 
from a stochastic simulation. All parameters are the 
same as in Fig. ||, the TM concentration is c = 10 /iM. 
The empty sites at the end of actin filaments are not 
counted as gaps, but they do contribute to the frac- 
tion of empty sites. The gap density in finite systems 
does not deviate significantly from infinite systems as 
long as it is higher than one gap per filament length. 



hilation at the ends) or raise (due to gap creation at the 
ends) the gap concentration. 



Time-dependent solution concentration 

In the previous sections we have assumed that the TM 
concentration c remains constant during the experiment. 
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FIG. 9: Average values of the gap density [na{t) 1 
dashed lines) and fraction of unoccupied binding sites 
(no(^)? solid lines) as a function of time for different 
actin concentrations ca- The simulation takes into ac- 
count the drop of TM concentration as part of it binds 
to actin. The results show that the plateau gap con- 
centration is not significantly influenced by this effect. 



This, however, was a simplification as the TM concentra- 
tion necessarily decreases when some of it gets bound to 
actin. Situations where TM is added gradually in course 
of the experiment are conceivable as well. 
Instead of Eq. ^| we now obtain 



fn{t) 

exp 



r n {c{t))P Q {t) =r„(c(t)) x 
r n (c(t')) ( f v(c(t"))dt" + 2d ) dt' 



(19) 



When a part of TM gets bound to actin, the solution 
concentration decreases according to (the expression for 
n (t) is analog to Eq. O) 



c(t) = c 



Co 



(l-no(i)) 



y(i -run 



3 / f n (t')dt' 
10 



(20) 



where Co denotes the initial TM concentration and ca the 
concentration of actin monomers forming the filaments. 
Equations (19) and ( |20|) uniquely determine the TM and 
gap concentration as a function of time and can be solved 
numerically. However, in most cases the difference to the 
solution with a constant concentration is not large. This 
is due to the fact that the gap concentration is deter- 
mined by the number of independent nucleation events 
and these mostly take place in the initial phase, when the 
solution concentration has not yet dropped significantly. 
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An example is shown in Fig. There the same curves 
as in Fig. f| are shown for the initial TM concentration of 
Co = 10 fiM and an actin concentration which can take up 
50% or 90% of all TM. Although the binding is somewhat 
slowed down, the plateau gap concentration stays practi- 
cally the same. In the second (gap annealing) stage the 
behavior is determined by the detachment rates, which 
are independent of the concentration. We therefore con- 
clude that the essential features of the system are cap- 
tured by the model with a constant TM concentration. 



DISCUSSION 

Table [n] shows the values of model parameters from 
the literature and the results of our calculation for these 
values. The equilibrium gap concentration n^ 9 and aver- 
age gap size g only depend on the binding constant K N , 
the cooperativity coefficient 7, and the TM solution con- 
centration c. For low concentrations (with a high gap 
concentration in equilibrium) the system reaches the final 
state very quickly. These situations are not the subject 
of our study and we therefore leave the fields in the table 
empty. We focus on cases with a high TM coverage of 
actin filaments. Interestingly, the gap concentration after 
the initial binding phase Uq only very weakly depends on 
the TM concentration and other model parameters, its 
value always being between 0.015 and 0.03 gaps per lat- 
tice site (between 2.8 and 5.6 gaps per micron). 

The gap concentration during the annealing phase de- 
pends almost entirely on the effective hopping rate r^ op , 
given by Eq. [l4|. The latter is of the order of magni- 
tude of the smaller among the two detachment rates for 
a molecule at the end of a block, and k_ . As sug- 



gested by Weigt et al. (1991 ), there is strong evidence for 
an asymmetry in the attachment rates on the plus (fc+) 
and the minus end (k+) of a block. We therefore esti- 
mate the hopping rate rh op between 0.1 s^ 1 and 0.6 s" 1 . 

In most cases listed in table |l| the equilibrium gap 
concentrations were very low. They are only relevant 
if the filaments are long enough to host at least a few 
gaps at this low concentration. In most experimental 
situations this will not be the case. We have shown that 
in a finite system the gap concentration follows the results 
for an infinite system as long as the gap concentration 
is at least a few per filament. We therefore estimate the 
relaxation time as the time in which the gap density drops 
below one per filament length. It scales as i Eq oc L 2 . For 
a filament of 1000 lattice sites (5.5 /mi) this gives a value 
between 1 and 6 hours. The relaxation time becomes at 
least 100 times shorter (40-220 s) if one considers a state 
with one gap per 100 lattice sites (roughly two gaps per 
micron) as equilibrated. 

The amount of data gathered in different experimen- 
tal studies allows us to make quite firm predictions for 
the vacancy density as a function of time. Our model- 
ing shows that in many relevant experimental situations, 
especially when using long filaments, the vacancy concen- 



tration is much higher than it would be in equilibrium. 
Actin filaments used in skeletal muscle sarcomeres, on 
the other hand, are short enough that the TM equili- 
brates within a minute. Therefore, care has to be taken 
when the calcium regulation of myosin activity is stud- 
ied on filaments assembled in vitro, as calcium regulation 
is by itself a strongly cooperative process and therefore 
vulnerable to gaps in TM filaments. We also show that 
with realistic parameters the effect of the initial solution 
concentration on the vacancy density is rather small (not 
more than a factor of 2). The only way to eliminate gaps 
is to give the system enough time to equilibrate (esti- 
mated as an hour on long filaments, though the gap con- 
centration already reaches quite low values after a few 
minutes), or to assemble the actin and the tropomyosin 
filaments simultaneously. 

Unfortunately, in none of the existing studies could the 
gap concentration be measured directly. But there could 
be indirect ways to test the predictions of our model. 
One possibility would be suddenly to decrease the TM 
concentration (by dilution) or binding affinity to actin 
(by an increase in the salt concentration) and study the 
detachment kinetics. As it is almost exclusively the TM 
molecules next to gaps that detach, the initial detach- 
ment rate could be a direct measure for the gap concen- 
tration. 

To conclude, we have shown that the dynamics of 
tropomyosin binding to actin polymers shows many 
stages. An initial fast binding phase is followed by a 
slow annealing process before it reaches the final equilib- 
rium. We could give analytical approximations for the 
time dependent gap density in all regimes, using elemen- 
tary nucleation-growth models for the first stage and and 
a mapping to a reaction-diffusion model in the second 
stage. 
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